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The effects of free-air carbon dioxide enrichment (FACE) and elevated soil and water temperature (warming) on the rice 
root-associated bacterial community were evaluated by clone library analysis of the 16S ribosomal RNA gene. Roots were 
sampled at the panicle initiation and ripening stages 41 and 92 days after transplanting (DAT), respectively. The relative 
abundances of the methanotrophs Methylosinus and Methylocystis were increased by warming and decreased by FACE at 92 
DAT, which indicated that microbial methane (CH 4 ) oxidation in rice roots may have been influenced by global warming. The 
relative abundance of Burkholderia kururiensis was increased by warming at 41 DAT and by FACE or warming at 92 DAT. 
The abundances of methanotrophs increased during rice growth, which was likely induced by an enhancement in the emission 
of CH 4 from the paddy fields, suggesting that CH 4 is one of the predominant factors affecting the structure of the microbial 
community in rice roots. Marked variations in the community structure were also observed during rice growth in other genera: 
Bradyrhizobium, Clostridium, and an unknown genus close to Epsilonproteobacteria were abundant at 92 DAT, whereas 
Achromobacter was abundant at 41 DAT. These results demonstrated that the community structures of rice root-associated 
bacteria were markedly affected by FACE, temperature, and the rice growth stage. 
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MICROBES AND 
ENVIRONMENTS 



The atmospheric carbon dioxide concentration ([C0 2 ]) 
was stable at 270 ppm for at least 1,000 years prior to the start 
of the industrial revolution. Since that time, [C0 2 ] has been 
rising and has reached nearly 400 ppm (7). The increase in 
[C0 2 ] is expected to enhance the growth and yield of C3 
crops, including rice (29). Previous studies reported that 
increases in [C0 2 ] quantitatively and qualitatively altered 
the release of labile sugars, organic acids, and amino acids 
from plant roots (2, 6), which may influence the activity of 
rhizospheric and root-associated microbes, including meth- 
anogenesis (18). Tokida et al. (32, 33) reported that the 
emission of CH 4 from paddy fields was significantly 
increased by [C0 2 ] and/or temperature elevations. CH 4 gen- 
erated in soil is diffused into rice roots, transported to the 
shoot via aerenchyma, and finally released from micropores 
in the leaf sheaths (25). The rhizosphere and rice roots were 
previously shown to be the main areas CH 4 oxidation in rice 
paddies (4) because the oxidation of CH 4 was inactive in 
flooded soils without oxygen derived from rice roots. 

Although the responses of rice plants to elevated atmos- 
pheric [C0 2 ] and/or temperature have been studied in detail, 
those of plant-associated microbes remain unknown (29, 33). 
The responses of plant-associated microbes to global climate 
changes may potentially be important because they play 
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major roles in the flow of carbon and nitrogen as the primary 
utilizers of plant-derived compounds in the rhizosphere. 

A technique to enrich bacterial cells obtained from plant 
tissues has provided a gateway to access plant-associated 
bacterial communities (14), and this has facilitated a deeper 
understanding of microbial community shifts caused by 
environmental factors (15, 27). In the present study, we 
investigated the community structure of rice root-associated 
bacteria in environments with elevated temperature (ET) and/ 
or [C0 2 ] using the bacterial cell enrichment method (14). Our 
results have provided an insight into carbon and nitrogen 
cycles in rice paddies under a changing climate. 

Materials and Methods 

Study site, [CO J enrichment, and soil and water warming 

The free-air carbon dioxide enrichment (FACE) and soil warming 
experiments were conducted in a rice paddy field at Tsukubamirai, 
Ibaraki, Japan (35°58'27"N, 139°59'32"E, 10 m above sea level), 
during the 201 1 growing season. The soil there is a Fluvisol, which 
is typical of alluvial areas. The bulk density was 0.87 Mg nr 3 , and 
total C and N were 21.4 mg g -1 and 1.97 mg g -1 , respectively. The 
cation exchange capacity was 202 umolc g -1 (12). The experimental 
site was established in 20 10, and control protocols for the FACE and 
warming treatments were described previously (23, 32). Briefly, 
four rice paddy fields were used as replicates, each with areas at 
ambient [C0 2 ] (AMBI) and also at enriched [COJ (FACE) with a 
target concentration of 200 ^mol moF 1 above AMBI. Each treat- 
ment area was a 240-m 2 octagon ("ring," hereafter). The FACE rings 
used emission tubes on all eight sides at a height of approximately 
30 cm above the canopy, and these released pure C0 2 from wind- 
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ward sides to maintain a stable concentration at the center of the 
rings (ambient + 200 mnol moL'). The AMBI and FACE rings were 
separated by at least 70 m (center to center), which was previously 
shown to be sufficient to prevent cross-contamination by CO-, from 
a FACE ring (13). 

Warming treatments were also conducted by a split-plot design in 
each ring with two levels of soil and water temperatures: normal 
temperature (NT) and ET with the target of 2°C above NT. Warming 
was achieved using heating wires placed on the soil surface between 
the rows, with the water temperature continuously measured by a 
PtlOO thermometer (Chino Co. Ltd., Tokyo, Japan). The water and 
plow layer (at a depth of 10 cm) temperatures were almost uniformly 
elevated. The ET plot was enclosed using corrugated PVC panels to 
prevent an exchange of the paddy water with the surrounding area. 

Rice cultivation and fertilization 

Rice (Oryza sativa L. cv. Koshihikari) was sown on 25 April 
2011 in seedling trays with 448 cells (Minora Pot 448, Minora 
Industrial Co., Ltd., Okayama, Japan). Three pre-germinated seeds 
were planted in each cell of the tray. After emergence, we raised 
seedlings in a puddled open field with a tunnel cloche or floating 
mulch for the first 2 weeks. On 25 and 26 May, seedlings at the 
five-leaf stage were transplanted into the rings by hand, with three 
seedlings per hill. Hills and rows were 15 and 30 cm apart, respec- 
tively, with a resultant density of 22.2 hills rrr 2 . Fertilizers were 
applied as a basal dressing. Nitrogen was supplied at 8 g N irr 2 (2 
and 6 g N irr 2 as urea and coated urea, respectively; 4 g of LP- 100 
and 2 g of LP-140; JCAM-Agri Co., Ltd., Tokyo, Japan). Phosphate 
and potassium were applied as a compound fertilizer (Sumitomo 
Chemical Co., Ltd., Tokyo, Japan) containing 4.4 g P irr 2 and 8.3 g 
K rrr 2 , respectively. The method for rice cultivation and fertilization 
was described previously (12). 

CH 4 emission measurements 

The emission of CH 4 was measured weekly or biweekly between 
7 June and 23 August using a closed chamber method, as described 
previously (16). Each chamber consisting of lower (60 cm H) and 
upper (60 cm H) sections was placed over 4 hills of rice plants with 
a basal area of 30 x 60 cm. The upper section of the chamber fit 
over the lower one and was supported by a water-filled groove 
surrounding the outer top lip of the lower section, thereby providing 
an airtight seal between the two sections and surrounding atmos- 
phere. Gas samples were collected from the chamber 0, 10, and 20 
min after placement of the chamber. The samples were injected into 
pre-evacuated 19-mL glass vials and brought back to the laboratory 
for analysis. The mixing ratio of CH 4 was determined by gas chro- 
matography equipped with a flame ionization detector (GC-14B; 
Shimazu, Kyoto, Japan). The emission of CH 4 was calculated based 
on an increase in the mixing ratio of the basal area of the chamber, 
chamber volume, and temperature inside the chamber. An analysis 
of variance (ANOVA) was conducted on the cumulative amount 
of CH 4 that was emitted during 41-90 days after transplanting 
(DAT) using a general linear procedure. [CO2], temperature, and 
[C0 2 ] x temperature were treated as fixed effects, while ring and 
ringx[C0 2 ] were treated as random effects. 

Rice sampling 

Three whole rice plants were carefully dug out by hand from each 
treatment plot on 5 July 2011 (41 d after transplanting [DAT]), 
which corresponded to the panicle initiation stage, and 25 August 
(92 DAT), which corresponded to the ripening stage (2 C0 2 treat- 
ments [AMBI or FACE] x 2 temperature treatments [ET or NT] x 4 
rings x 2 growth stages). A block of plow-layer soil, 30 x 15 x 30 cm 
(length x width x depth), was also taken with the plants. Whole rice 
plants were then immediately transported to the laboratory and 
washed with tap water to remove adhering soil particles. The roots 
were separated from the shoot and stored at -80°C until they were 
used for analysis. 



DNA preparation and clone libraiy construction 

The root samples were manually ground into a fine powder in 
liquid nitrogen using a mortar and pestle. Three ground root 
samples, collected from the same ring of the same treatment, were 
composited and then homogenized again in a blender. The com- 
posited samples were subjected to DNA extraction by the bacterial 
cell enrichment method (14), PCR amplification, and clone library 
analysis targeting the 16S rRNA gene (14). Briefly, bacterial DNA 
was extracted from the composited roots using the bacterial cell 
enrichment method (14). The PCR clone libraries for 16S rRNA 
genes were constructed as follows: 10 ng total bacterial DNA was 
used as a template in a final reaction volume of 50 p.L, including 0.1 
uM of each primer and 2 U of Ex Taq DNA polymerase (Takara 
Bio, Otsu, Japan) with the universal primers 27F (5'-AGAGTTTGAT 
CMTGGCTCAG-3') and 1525R (5'-AAGGAGGTGWTCCARCC 
-3') (17). The cycling conditions were an initial denaturation step of 
2 min at 94°C; 25 cycles of 30 s at 94°C, 30 s at 55°C, and 1.5 min 
at 72°C; and a final extension step of 8 min at 72°C. PCR products 
of the predicted size (1500 bp) were purified using a Wizard SV Gel 
and PCR Clean-Up System (Promega Japan, Tokyo, Japan). Four 
composited samples, collected from four rings within the same 
treatment, were combined, then ligated into the pGEM-T Easy 
plasmid vector (Promega Japan) at 25°C for 1 h. A partial sequence 
of the 16S rRNA gene was determined by Takara Bio Inc. (Otsu, 
Japan) using the 27F forward primer as a sequencing primer. 

Bioinformatics analysis 

Low-quality and contaminated reads were removed using 
OrientationChecker (1). A partial 16S rRNA gene sequence (corre- 
sponding to bases 109 to 684 of the Escherichia coli 16S rRNA 
gene) was used for sequence analyses. Chimeric sequences were 
detected using MALLARD (1). The remaining sequences were 
aligned using CLUSTAL W (31). Based on the alignment, a distance 
matrix was constructed using the DNADIST program from PHYLIP 
ver. 3.66 (http://evolution.genetics.washington.edu/phylip.html) with 
the default parameters. The resulting matrix was clustered using 
Mothur (28) to generate diversity indexes with a threshold of a 97% 
sequence identity. Library coverage was calculated with the nonpara- 
metric estimator C (9). UniFrac (19) was applied with the abundance- 
weighted option to examine similarities among clone libraries. 

The phylogenetic composition of the library sequences was eval- 
uated using the RDPmultiple classifier (34), with confidence levels 
of 80%. Sequences assigned to Burkholderia and Bradyrhizobium 
were extracted separately and aligned with reference sequences 
using CLUSTAL W (31). Neighbor-joining trees were constructed 
using MEGA version 5.1 (30), and 1,000 bootstrap replicates were 
used to generate a consensus tree. 

Nucleotide sequence accession numbers 

The nucleotide sequences of 16S rRNA genes in the clone 
libraries have been deposited in DDBJ under the accession numbers 
shown in Table 1 . 

Results and Discussion 

Overview of bacterial community structures 

The statistics of the clone libraries are summarized in 



Table 1. DDBJ accession numbers of 16S rRNA gene sequences 



Sample 


Panicle initiation stage 
(41 DAT) 


Ripening stage 
(92 DAT) 


AMBI-NT 
AMBI-ET 
FACE-NT 
FACE-ET 


AB836880-AB837055 
AB837056-AB837234 
AB837235-AB837406 
AB837407-AB837584 


AB837585-AB837749 
AB837750-AB837919 
AB837920-AB838080 
AB838081-AB838242 



DAT, days after transplanting; AMBI, ambient C0 2 ; FACE, free-air 
C0 2 enrichment; NT, normal soil and water temperature; ET, elevated 
soil and water temperature. 
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Table 2. Statistical characteristics of 16S rRNA gene clone libraries derived from rice roots 
Rice growth stage Panicle initiation stage (41 DAT) Ripening stage (92 DAT) 



C0 2 AMBI FACE AMBI FACE 



Temperature 


NT 


ET 


NT 


ET 


NT 


ET 


NT 


ET 


Statistics 


















No. of sequences 


176 


179 


172 


178 


165 


170 


161 


162 


No. of OTUs (>97% identity) 


28 


22 


27 


28 


52 


45 


48 


32 


No. of singletons 


20 


15 


19 


20 


35 


28 


30 


17 


Library coverage (%) a 


88.6 


91.6 


89.0 


88.8 


78.8 


83.5 


81.4 


89.5 


Diversity indexes 


















Chaol 


91.3 


57.0 


69.8 


75.5 


118.1 


82.8 


96.3 


66.0 


Shannon 


2.2 


1.5 


2.1 


1.9 


3.3 


3.0 


3.1 


2.7 



DAT, days after transplanting; AMBI, ambient C0 2 ; FACE, free-air CO2 enrichment; NT, normal soil and water temperature; ET, elevated soil and 
water temperature; OTU, operational taxonomic unit. 

a Coverage calculated as C x = 1 - (n x IN), where n x is the number of singletons that are encountered only once in a library and N is the total number 
of clones. 
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Fig. 1. Principal coordinate analysis of the 16S rRNA gene libraries 
of bacterial communities in rice roots under normal and elevated [CO2] 
and temperature conditions. The ordination was constructed using 
UniFrac distances weighted by the relative abundances. Principal 
component 1 (PC 1 ) and principal component 2 (PC2) are plotted on the 
x- and>'-axes, respectively. The percentage of variation explained by the 
plotted principal coordinates is indicated on the axes. Samples were 
collected at the panicle initiation stage (•, 41 DAT) and ripening stage 
(■, 92 DAT). 



Table 2. In all treatments, the number of operational taxo- 
nomic units (OTUs) and the Chaol and Shannon diversity 
indexes were greater at 92 DAT than at 41 DAT. These 
indexes were decreased in the samples at 92 DAT due to the 
elevation in [C0 2 ] or soil and water temperature. No clear 
trend was observed in the samples at 41 DAT. 

Principal coordinate analysis was performed using all 
sequence data (Fig. 1) in order to obtain an overview of 
bacterial community shifts caused by the rice growth stage 
and elevation in [C0 2 ] and temperature. Samples were clearly 
separated along the first principal component (PCI) axis 
(64.86%) according to the rice growth stage, which indicated 
that bacterial community structures markedly changed as the 
host plant grew. Two tight clusters were formed according to 



the temperature condition in samples at 4 1 DAT, suggesting 
that community structures were more sensitive to the 
temperature change than to that of [C0 2 ] at 41 DAT. 
Community shifts in samples at 92 DAT were more com- 
plicated. The degree of the community shift from the control 
(AMBI-NT) was smaller in the simultaneous treatment 
(FACE-ET) than in the other treatments (FACE-NT and 
AMBI-ET). Furthermore, FACE-NT and AMBI-ET were 
clustered close to each other at 92 DAT. 

Phylogenetic composition 

In all treatments, the abundances of Methylosinus and 
Methylocystis were markedly higher at 92 DAT (5.0-15.3%) 
than at 41 DAT (0.0-1.1%) (Table 3). At the beginning of the 
rice growing period, the amount of CH 4 emitted was very low 
(Table 4). As the season progressed, it steadily increased to 
approximately 14-18 mg C-CH 4 nr 2 h" 1 at 48 DAT (July 12), 
and a high emission level of CH 4 was maintained until 83 
DAT (August 16). These results suggest that rice roots were 
exposed to a large amount of CH 4 between the first (41 DAT) 
and second sampling (92 DAT), and this may have caused the 
increase observed in the relative abundances of Methylosinus 
and Methylocystis in the rice roots at 92 DAT. The cumula- 
tive emission of CH 4 during 41-90 DAT was the highest in 
FACE-ET (Table 4). However, the relative abundances of 
Methylosinus and Methylocystis in FACE-ET were inter- 
mediate among the four treatments (Table 3), which sug- 
gested that factors other than CH 4 also affected the relative 
abundances of methanotrophs. At 92 DAT, the relative 
abundances of Methylosinus and Methylocystis were 
increased by the elevation in temperature in both AMBI and 
FACE plots, but were decreased by that in [C0 2 ]. A previous 
FACE experiment conducted in Japan showed that nitrogen 
concentrations in rice plants were decreased by elevations in 
[C0 2 ] (29), which was at least partially attributed to a dilution 
effect due to the greater production of dry matter. Many 
studies have suggested the stimulatory effect of nitrogen 
on CH 4 oxidation in rice paddies (3). Such a change in 
the nitrogen condition may affect the activities of rice root- 
associated methanotrophs, leading to decreases in the relative 
abundance by elevations in [C0 2 ]. Tokida et al. (32) previ- 
ously reported that the emission of CH 4 from paddy fields 
was significantly increased by elevations in [C0 2 ], and this 
effect was considered to be mainly derived from an increase 
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Table 3. Phylogenetic compositions of 16S rRNA gene clone libraries derived from rice roots (values are relative abundances as a percentage of 
total bacteria) 



Rice growth sttigc 


Panicle initiation stage (41 DAT) 




Ripening 


stase (92 DAT) 




co 2 


AMBI 


FACE 


AMBI 


FACE 




1 emperature 


NT ET 


NT ET 


NT 


ET 


VT 
IN 1 


FT 


Proteobacteria 














Alphaproteobacteria 


10.8 3.4 


5.2 5.6 


17.6 


26.5 


15.5 


12.3 


Methylosinus (A) 


0.0 0.0 


0.0 0.0 


6.1 


11.2 


2.5 


4.3 


Methylocystis (B) 


0.0 0.6 


0.0 1.1 


3.0 


4.1 


2.5 


2.5 


(A)+(B) 


0.0 0.6 


0.0 1.1 


9.1 


15.3 


5.0 


6.8 


Bradyrhizobium 


0.0 0.6 


1.2 1.1 


5.5 


4.1 


5.6 


4.9 


Rhizobium 


10.2 2.2 


4.1 1.7 


1.8 


4.7 


1.9 


0.0 


Betaproteobacteria 


77.8 87.2 


77.9 84.3 


10.9 


42.4 


47.8 


35.8 


Burkholderia 


44.9 72.1 


34.9 64.0 


8.5 


37.6 


42.9 


35.2 


Ralstonia 


27.8 13.4 


39.5 15.7 


2.4 


3.5 


4.3 


0.6 


Achromobacter 


4.0 1.7 


2.3 3.9 


0.0 


0.0 


0.0 


0.0 


Gammaproteobacteria 


2.3 1.1 


3.5 2.8 


5.5 


2.4 


2.5 


6.2 


Aquicella 


0.6 0.6 


1.2 1.7 


5.5 


1.8 


2.5 


5.6 


Deltaproteobacteria 


0.0 0.0 


0.0 0.0 


1.2 


0.6 


0.0 


0.0 


Firmicutes 














Bacilli 


1.1 0.0 


0.6 0.6 


0.0 


0.0 


0.0 


0.0 


Clostridia 


1.7 0.0 


1.2 1.1 


9.7 


5.9 


11.2 


13.6 


Clostridium III 


0.0 0.0 


0.0 0.0 


7.3 


3.5 


3.1 


9.9 


Clostridium sensu stricto 


0.0 0.0 


0.6 0.6 


0.0 


0.0 


3.1 


1.2 


Others 


1.1 2.2 


0.6 0.6 


1.8 


1.8 


2.5 


0.6 


Unclassified 


5.1 6.1 


11.0 5.1 


53.3 


20.6 


20.5 


31.5 


Unclassified OTU164 


0.0 0.0 


0.0 0.0 


33.9 


1.2 


2.5 


18.5 


DAT, days after transplanting; AMBI, ambient CO2; FACE, free- 


air CO2 enrichment; NT, normal soil and water temperature; ET, elevated soil and 


water temperature. 














Table 


4. Seasonal and total CH4 emission from each rice paddy treatment 








AMBI-NT 


FACE-NT 


AMBI-ET 


FACE-ET 


13 DAT (7 June) 


0.64 ±0.16 


0.81 ±0.16 


0.64 ±0.19 


0.51 ±0.42 


27 DAT (21 June) 


1.88 ±0.81 


2.68 ±0.69 


2.85 ±1.68 


3.47 ± 1.87 


41 DAT (5 July) 


7.85 ± 1.62 


10.87 ±2.98 


9.96 ±2.82 


11.98 ±2.67 


48 DAT (12 July) 


14.00 ±2.23 


15.89 ±3.38 


17.26 ±2.78 


18.35 ±3.14 


62 DAT (26 July) 


14.71 ±3.97 


14.39 ±3.61 


18.06 ±4.08 


18.69 ±2.40 


69 DAT (2 August) 


15.01 ±7.80 


16.25 ±4.15 


15.42 ±3.77 


18.37 ±5.09 


76 DAT (9 August) 


30.36 ±8.50 


22.36 ± 6.44 


27.50 ±7.63 


27.16 ±6.63 


83 DAT (16 August) 


22.98 ±4.55 


19.96 ±5.44 


26.72 ±9. 12 


31.83 ±8.1 


18 


90 DAT (23 August) 


11.71 ±5.70 


14.96 ± 10.29 


13.30 ±4.62 


14.33 ±9.67 


Total emission (41-90 DAT) 


20.36 ± 5.04 


19.64 ±4.69 


22.22 ± 2.99 


24.13 ±3.33 


ANOVA results 














[C0 2 ] 




not significant (P=0.67) 








Temperature 




P=0.05 










[CO2] xTemperature 




not significant (P=0.34) 








Data are the mean ± SD of four replicated plots (except n = 3 for AMB1-ET at 62 DAT and FACE-ET at 83 DAT). 






Seasonal CH 4 emissions are shown in mg 


C-CH4 m~ 2 fr 1 and total CH 4 emission in g C-CH4 m~ 2 . 










Total emission shows the cumulative emission of CH 4 during 41- 


-90 DAT, which almost corresponded to the two microbial sampling dates. 


DAT, days after transplanting; AMB1, ambient CO2; FACE, free-air CO2 enrichment; NT, normal soil and water temperature; ET, elevated soil and 



water temperature. 



in rhizodeposition. However, our results suggest that one 
reason for the increase in CH 4 emission with elevations in 
[C0 2 ] may have been a decline in the abundance of methano- 
trophs associated with rice roots. An enhancement in CH 4 
oxidation activity in rice roots is vital for breaking the 
positive feedback loop of CH 4 emission that will occur with 
increases in atmospheric [C0 2 ] in the future. However, we 
did not observe a clear increase in the emission of CH 4 by 
elevations in [C0 2 ] under the NT condition (Table 4). 
Therefore, the effects of [C0 2 ] and temperature elevation 
on the oxidation activity of CH 4 in paddy fields need to be 
studied in more detail. 

Members of the genus Bradyrhizobium are important 



nitrogen-fixing bacteria in rice roots (5). In all treatments, the 
relative abundance of Bradyrhizobium was greater at 92 DAT 
(4.1-5.6%) than at 41 DAT (0.0-1.2%) (Table 3). However, 
no apparent effects of the elevation in [C0 2 ] and temperature 
were observed at 41 or 92 DAT. To perform a detailed phylo- 
genetic analysis of rice root-associated Bradyrhizobium, 16S 
rRNA reads assigned to Bradyrhizobium were extracted from 
clone libraries and a phylogenetic tree was constructed with 
the other members of bradirhizobia (Fig. 2A). Rice root- 
associated bradyrhizobia were clustered into two groups 
that were phylogenetically close to Bradyrhizobium sp. 
ORS278 (bradyrhizobial cluster I) (8, 26) and Bradyrhizobium 
jicamae (bradyrhizobial cluster II). Bradyrhizobial cluster 



188 



OKUBO et al. 



S 3 



a 1 



bradyrhizobial cluster I 

19 sequences (Fig. 2B) 

Bradyrhizobium sp. ORS278 (AF239255) 

58,2FNT1_B05 (FACE NT, Ripening, 92 DAT) 



Bradyrhizobium oiigotrophicumS5B (D78366) 

84 , Bradyrhizobium cytisiCT AW/1 1 (EU561 065) 



1 Bradyrhizobium betae LMG21987 (AY372184) 
Bradyrhizobium iriomotenseEK05 (AB300992) 
Bradyrhizobiumsp. CCBAU 23303 (HQ231463) 
'' Bradyrhizobium denitrificans LMG8443 (X66025) 

Bradyrhizobium japonicum USDA1 1 0 (BA000040) 

46 I Bradyrhizobium japonicum LMG6138 (X66024) 

' Bradyrhizobiumsp. BTA-1 (AJ558025) 

Bradyrhizobium liaoningenseLUGI&Zid (AJ250813) 
Bradyrhizobium yuanmingense CCBAU 1 1 07 1 (AF1 9381 8) 

- Bradyrhizobium jicamaePAC68 (AY624134) 

bradyrhizobial cluster II 

18 sequences (Fig. 2C) 



- Bradyrhizobium /ab/abi CCBAU23086 (GU433448) 
i Bradyrhizobium e/kaniiUSDA76 (U35000) 



87 ' Bradyrhizobium pachyrhiziP AC4& (AY624 135) 
Bosea thiooxidansDSM9653 (AJ250796) 



bradyrhizobial 
cluster I 
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lllii ■■■hll 
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Fig. 2. Phylogenetic position and distribution of rice root-associated Bradyrhizobium. (A) Phylogenetic tree of rice root-associated bradyrhizobia 
extracted from the clone libraries. (B, C) Distributions of reads belonging to the two large groups in panel (A). Accession numbers are shown 
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I was only observed at 92 DAT (Fig. 2B), whereas bradyrhi- 
zobial cluster II was observed at both 41 and 92 DAT (Fig. 
2C). In a previous study, Bradyrhizobium sp. ORS278 was 
reported to colonize the surface and intercellular space of rice 
roots and also fix nitrogen (5). Our clone library analysis 



suggested that members of bradyrhizobial cluster I may be 
representative nitrogen-fixing bacteria in the rice root at 
the ripening stage. The relative abundance of bradyrhizobial 
cluster I (Fig. 2B) was strongly correlated with those of 
Methylosinus and Methylocystis (Pearson's correlation coef- 
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ficient, r=0.83, P=0.011; r=0.92, P=0.001, respectively). 
One explanation for this correlation was the possible meta- 
bolic coupling by which one-carbon compounds oxidized by 
methanotrophs such as methanol may be partially consumed 
by these Bradyrhizobium species. Bradyrhizobium sp. 
ORS278 (bradyrhizobial cluster I) has a gene for methanol 
oxidation (BRAD05483-BRAD05487) on the genome, and 
was able to oxidize methanol (Seki et al. unpublished result). 

Burkholderia was the most dominant genus in all samples, 
except AMBI-NT at 92 DAT, (Table 3). Members of this 
genus differ in terms of their effects on rice plants by 
exhibiting pathogenic {e.g., Burkholderia glumae) (11) or 
symbiotic {e.g., Burkholderia kururiensis) interactions (21). 
A phylogenetic tree was constructed using 16S rRNA reads 
assigned to rice root-associated Burkholderia in the clone 
libraries (Fig. 3A). Most of the sequences clustered into 
one group that was phylogenetically close to B. kururiensis 
(Figs. 3A, B). This species was reported to colonize rice roots 
and significantly enhance rice growth by fixing nitrogen and 
producing the phytohormone auxin (21). Previous rice FACE 
experiments showed that the nitrogen concentration of rice 
was decreased by elevations in [C0 2 ] (29), which suggested 
that nitrogen availability is a limiting factor in an elevated 
[C0 2 ] environment. In the present study, the relative abun- 
dance of Burkholderia was still high at 92 DAT in under the 
elevated [C0 2 ] condition (Table 3). These results suggest that 
when a nitrogen deficiency occurs in rice, the high relative 
abundance of Burkholderia is maintained for a longer period 
in order to support growth. The high relative abundance of 
Burkholderia was also maintained at 92 DAT in the ET 
condition (Table 3). Root-associated Burkholderia may also 
play roles to support the growth of rice stimulated by eleva- 
tions in temperature (33). 

A high relative abundance of Clostridia was observed in 
all treatments at 92 DAT (5.9-13.6%) (Table 3), suggesting 
the presence of an anaerobic environment in rice roots at 
the ripening stage. At 92 DAT, Clostridium cluster III was 
the most abundant genus in Clostridia, the members of 
which produce cellulosomes (highly active cellulolytic and 
xylanolytic complexes) (24). Plant biomass decomposed by 
Clostridium cluster III and other microbes may be a substrate 
for rice root-associated methanogenic archaea (15, 18). 

One large OTU (OTU164), which was an unclassified 
bacterium showing 89% similarity to Nitratiruptor tergarcus, 
was abundant in AMBI-NT (33.9%) and FACE-ET 
(18.5%) at 92 DAT (Table 3). Nitratiruptor tergarcus is 
an Epsilonproteobacteria and a nitrate-reducing chemo- 
lithoautotroph (22). This OTU member is one of the main 
reasons why AMBI-NT and FACE-ET at 92 DAT were 
clustered with each other in the principal coordinate analysis 
(Fig. 1). The relative abundance of OTU164 was markedly 
decreased by the elevation in temperature (1.2% vs. 33.9% in 
AMBI-NT) or [C0 2 ] (2.5%), but was only slightly decreased 
by the simultaneous elevation in temperature and [C0 2 ] 
(18.5%) (Table 3). These results suggested that the effects of 
temperature and [C0 2 ] elevations interact with each other and 
changes in the community structure are compensated. 

Achromobacter, a genus of sulfur-oxidizing bacteria in 
rice paddy fields (10), was only observed at 41 DAT (Table 
3). Sulfur-oxidizing bacteria are able to oxidize reduced 



sulfur compounds (H 2 S, thiosulfate, and sulfite). Sulfate is 
the main sulfur compound that rice roots can take up (20). 
Achromobacter may play an important role in the growth of 
rice at the panicle initiation stage. 

Conclusion 

[C0 2 ] enrichment (FACE), elevated soil & water tempera- 
ture (warming), and rice growth stages markedly affected 
the microbial communities of rice root-associated bacteria 
including Methylosinus sp., Methylocystis sp., Burkholderia 
kururiensis, Bradyrhizobium sp., Clostridium sp., and an 
unknown genus (OTU164) close to Epsilonproteobacteria. 
Most of these bacteria play important roles in the metabolism 
of C and N in the environment through, for example, nitrogen 
fixation and methane oxidation. The results of the present 
study will contribute to improving our understanding of 
microbe-mediated CN dynamics in paddy rice fields under a 
changing climate. 

Acknowledgements 

This work was supported in part by grants from the Japan Society 
for the Promotion of Science (JSPS) [grant number H2247071, 
H3247071, and H260430] to TO; the Ministry of Agriculture, 
Forestry and Fisheries of Japan (Genomics for Agricultural 
Innovation PMI-0002, and BRAIN); and the Ministry of Education, 
Culture, Sports, Science and Technology of Japan (Grants-in-Aid 
for Scientific Research (A) 23248052 to KM). 

References 

1. Ashelford, K.E., NA. Chuzhanova, J.C. Fry, A.J. Jones, and A.J. 
Weightman. 2006. New screening software shows that most recent 
large 16S rRNA gene clone libraries contain chimeras. Appl. Environ. 
Microbiol. 72:5734-5741. 

2. Bhattacharyya, P., K.S. Roy, S. Neogi, M.C. Manna, T.K. Adhya, 
K.S. Rao, and A.K. Nayak. 2013. Influence of elevated carbon dioxide 
and temperature on belowground carbon allocation and enzyme 
activities in tropical flooded soil planted with rice. Environ. Monit. 
Assess. 185:8659-8671. 

3. Bodelier, P.L., and H.J. Laanbroek. 2004. Nitrogen as a regulatory 
factor of methane oxidation in soils and sediments. FEMS Microbiol. 
Ecol. 47:265-277. 

4. Bosse, U., and P. Frenzel. 1997. Activity and distribution of 
methane-oxidizing bacteria in flooded rice soil microcosms and in 
rice plants {Oryza sativa). Appl. Environ. Microbiol. 63: 1 199-1207. 

5. Chaintreuil, C, E. Giraud, Y. Prin, J. Lorquin, A. Ba, M. Gillis, P. de 
Lajudie, and B. Dreyfus. 2000. Photosynthetic bradyrhizobia are 
natural endophytes of the African wild rice Oryza breviligulata. 
Appl. Environ. Microbiol. 66:5437-5447. 

6. Drigo, B., GA. Kowalchuk, and J.A. van Veen. 2008. Climate change 
goes underground: effects of elevated atmospheric C0 2 on microbial 
community structure and activities in the rhizosphere. Biol. Fertil. 
Soils 44:667-679. 

7. Forster, P., V. Ramaswamy, P. Artaxo, et al. 2007. Changes in 
atmospheric constituents and in radiative forcing, p. 129-234. In S. 
Solomon, D. Qin, M. Manning, Z. Chen, M. Marquis, K.B. Averyt, M. 
Tignor, and H.L. Miller, Climate Change 2007: The Physical Science 
Basis, Contribution of Working Group I to the Fourth Assessment 
Report of the Intergovernmental Panel on Climate Change. Cambridge 
University Press, Cambridge, UK. 

8. Giraud, E., L. Moulin, D. Vallenet, et al. 2007. Legumes symbioses: 
absence of Nod genes in photosynthetic bradyrhizobia. Science 
316:1307-1312. 

9. Good, I.J. 1953. The population frequencies of species and the estima- 
tion of population parameters. Biometrika 40:237-264. 

10. Graff, A., and S. Stubner. 2003. Isolation and molecular characteriza- 
tion of thiosulfate-oxidizing bacteria from an Italian rice field soil. 
Syst. Appl. Microbiol. 26:445^152. 



190 



OKUBO et al. 



11. Ham, J.H., R.A. Melanson, and M.C. Rush. 2011. Burkholderia 
glumae: Next major pathogen of rice? Mol. Plant Pathol. 12:329-339. 

12. Hasegawa, T., H. Sakai, T. Tokida, et al. 2013. Rice cultivar responses 
to elevated C0 2 at two free-air C0 2 enrichment (FACE) sites in Japan. 
Funct. Plant Biol. 40:148-159. 

13. Heim, A., N. Moser, H. Blum, and M.W.I. Schmidt. 2009. How far do 
experimentally elevated C0 2 levels reach into the surrounding? An 
example using the 13 C label of soil organic matter as an archive. 
Global Change Biol. 15:1598-1602. 

14. Ikeda, S., T. Kaneko, T. Okubo, et al. 2009. Development of a 
bacterial cell enrichment method and its application to the community 
analysis in soybean stems. Microbial. Ecol. 58:703-714. 

15. Ikeda, S., K. Sasaki, T. Okubo, et al. 2014. Low nitrogen fertilization 
adapts rice root microbiome to low nutrient environment by changing 
biogeochemical functions. Microbes Environ. 29:50-59. 

16. Inubushi, K., W. Cheng, S. Aonuma, M.M. Hoque, K. Kobayashi, 
S. Miura, H.Y. Kim, and M. Okada. 2003. Effect of free-air C0 2 
enrichment (FACE) on CH 4 emission from a rice paddy field. Global 
Change Biol. 9:1458-1464. 

17. Lane, D.J. 1991. 16S/23S rRNA sequencing, p. 115-175. In E. 
Stackebrandt and M. Goodfellow (ed.), Nucleic Acid Techniques in 
Bacterial Systematics. Wiley, Chichester, UK. 

18. Liu, G.C., T. Tokida, T. Matsunami, H. Nakamura, M. Okada, R. 
Sameshima, T. Hasegawa, and S. Sugiyama. 2012. Microbial commu- 
nity composition controls the effects of climate change on methane 
emission from rice paddies. Environ. Microbiol. Rep. 4:648-654. 

19. Lozupone, C, and R. Knight. 2005. UniFrac: a new phylogenetic 
method for comparing microbial communities. Appl. Environ. 
Microbiol. 71:8228-8835. 

20. Lunde, C, A. Zygadlo, H.T. Simonsen, P.L. Nielsen, A. Blennow, 
and A. Haldrup. 2008. Sulfur starvation in rice: the effect on photo- 
synthesis, carbohydrate metabolism, and oxidative stress protective 
pathways. Physiol. Plant 134:508-521. 

21. Mattos, K.A., V.L. Padua, A. Romeiro, et al. 2008. Endophytic 
colonization of rice (Oryza sativa L.) by the diazotrophic bacterium 
Burkholderia kururiensis and its ability to enhance plant growth. Ann. 
Acad. Bras. Cienc. 80:477-493. 

22. Nakagawa, S., K. Takai, F. Inagaki, K. Horikoshi, and Y. Sako. 
2005. Nitratiruptor tergarcus gen. nov., sp. nov. and Nitratifractor 
salsuginis gen. nov., sp. nov., nitrate-reducing chemolithoautotrophs 
of the epsilon-Proteobacteria isolated from a deep-sea hydrothermal 
system in the Mid-Okinawa Trough. Int. J. Syst. Evol. Microbiol. 
55:925-933. 

23. Nakamura, H, T. Tokida, M. Yoshimoto, H. Sakai, M. Fukuoka, and 
T. Hasegawa. 2012. Performance of the enlarged rice-FACE system 
using pure C0 2 installed in Tsukuba, Japan. J. Agric. Meteorol. 
68:15-23. 



24. Nishiyama, T., A. Ueki, N. Kaku, and K. Ueki. 2009. Clostridium 
sufflavum sp. nov., isolated from a mefhanogenic reactor treating cattle 
waste. Int. J. Syst. Evol. Microbiol. 59:981-986. 

25. Nouchi, I., S. Mariko, and K. Aoki. 1990. Mechanism of methane 
transport from the rhizosphere to the atmosphere through rice plants. 
Plant Physiol. 94:59-66. 

26. Okubo, T., S. Fukushima, M. Itakura, et al. 2013. Genome analysis 
suggests that the soil oligotrophic bacterium Agromonas oligotrophica 
(Bradyrhizobium oligotrophicum) is a nitrogen-fixing symbiont of 
Aeschvnomene indica. Appl. Environ. Microbiol. 79:2542-2551. 

27. Sasaki, K., S. Ikeda, T. Ohkubo, C. Kisara, T. Sato, and K. 
Minamisawa. 2013. Effects of plant genotype and nitrogen level on 
bacterial communities in rice shoots and roots. Microbes Environ. 
28:391-395. 

28. Schloss, P.D., S.L. Westcott, T. Ryabin, et al. 2009. Introducing 
Mothur: open-source, platform-independent, community-supported 
software for describing and comparing microbial communities. Appl. 
Environ. Microbiol. 75:7537-7541. 

29. Shimono, H, M. Okada, Y. Yamakawa, H. Nakamura, K. Kobayashi, 
and T. Hasegawa. 2009. Genotypic variation in rice yield enhance- 
ment by elevated C0 2 relates to growth before heading and not to 
maturity group. J. Exp. Bot. 60:523-532. 

30. Tamura, K, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. 
Kumar. 2011. MEGA5: molecular evolutionary genetics analysis 
using maximum likelihood, evolutionary distance, and maximum 
parsimony methods. Mol. Biol. Evol. 28:2731-2739. 

31. Thompson, J.D., D.G. Higgins, and T.J. Gibson. 1994. CLUSTAL W: 
improving the sensitivity of progressive multiple sequence alignment 
through sequence weighting, position-specific gap penalties and 
weight matrix choice. Nucleic Acids Res. 22:4673-4680. 

32. Tokida, T., T. Fumoto, W. Cheng, et al. 2010. Effects of free-air C0 2 
enrichment (FACE) and soil warming on CH 4 emission from a rice 
paddy field: impact assessment and stoichiometric evaluation. 
Biogeosciences 7:2639-2653. 

33. Tokida, T., M. Adachi, W. Cheng, et al. 201 1 . Methane and soil C0 2 
production from current-season photosynthates in a rice paddy 
exposed to elevated COi concentration and soil temperature. Global 
Change Biol. 17:3327-3337. 

34. Wang, Q., G.M. Garrity, J.M. Tiedje, and J.R. Cole. 2007. Naive 
Bayesian classifier for rapid assignment of rRNA sequences into the 
new bacterial taxonomy. Appl. Environ. Microbiol. 73:5261-5267. 



